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ABSTRACT 

The CO(l-O) and (2-1) emission of the circumstellar envelope of the Asymptotic Giant Branch (AGB) star EP Aqr has been observed 
in 2003 using the IRAM Plateau-de-Bure Interferometer and in 2004 using the IRAM 30-m telescope at Pico Veleta. The line profiles 
reveal the presence of two distinct components centered on the star velocity, a broad component extending up to ~10 kms“* and 
a narrow component indicating an expansion velocity of only ~2 kms^'. An early analysis of these data was performed under the 
assumption of isotropic winds. The present study revisits this interpretation by assuming instead a bipolar outflow nearly aligned with 
the line of sight. A satisfactory description of the observed flux densities is obtained with a radial expansion velocity increasing from 
~2 kms“* at the equator to ~10 kms^' near the poles. The mass loss rate is ~1.8 10“^ Moyr“*. The angular aperture of the bipolar 
outflow is ~45° with respect to the star axis, which makes an angle of ~13° with the line of sight. A detailed study of the CO(l-O) to 
CO(2-l) flux ratio reveals a significant dependence of the temperature on the star latitude, smaller and steeper at the poles than at the 
equator at large distances from the star (> 2" = 1.0 x 10“^ pc). Under the hypothesis of radial expansion of the gas and of rotation 
invariance about the star axis, the effective density has been evaluated in space as a function of star coordinates (longitude, latitude 
and distance from the star). Evidence is found for an enhancement of the effective density in the northern hemisphere of the star at 
angular distances in excess of ~3" and covering the whole longitudinal range. The peak velocity of the narrow component is observed 
to vary slightly with position on the sky, a variation consistent with the model and understood as the effect of the inclination of the 
star axis with respect to the line of sight. This variation is inconsistent with the assumption of a spherical wind and strengthens our 
interpretation in terms of an axisymmetric outflow. While the phenomenological model presented here reproduces well the general 
features of the observations, not only qualitatively hut also quantitatively, significant differences are also revealed, which would require 
a better spatial resolution to be properly described and understood. 

Key words. Stars: AGB and post-AGB - (Stars:) circumstellar matter - Stars; individual: EP Aqr - Stars; mass-loss - radio lines: 
stars. 


1. Introduction 

EP Aqr is one of the nearest mass losing Asymptotic Giant 
Branch (AGB) stars (d-WA pc, van Leeuwen 2007), and as such 
one of the best characterized objects of its class. The absence 
of Technetium in the spectrum (Lebzelter & Hron 1999), and 
the low value of the abundance ratio (~ 10, Cami et al. 

2000) indicate that it is still at the beginning of its evolution on 
the AGB. Yet Herschel has imaged a large scale (2' x4') circum¬ 
stellar shell at 70 pm (Cox et al. 2012) testifying for a relatively 
long duration of the mass loss (>10^ years) and for interaction 
with the surrounding interstellar medium. The elongation of the 
far-infrared image is in a direction opposite to the EP Aqr space 
motion. Prom Hi observations at 21 cm, Le Bertre & Gerard 
(2004) estimate a duration of ~ 1.6x10^ years for the present 
episode of mass loss. 

Combined observations have been obtained with the IRAM 
30-m telescope and the Plateau-de-Bure Interferometer in 
^^CO(l-O) and ^^CO(2-l) in 2003 and 2004, improving over 
earlier observations (Knapp et al. 1998, Nakashima 2006). The 


* Based on observations carried out with the IRAM Plateau-de-Bure 
Interferometer and the IRAM 30-m telescope. IRAM is supported by 
INSU/CNRS (France), MPG (Germany) and IGN (Spain). 


spectra, which reveal a wind composed of two main compo¬ 
nents with expansion velocities ~2 and ~10 kms"', were ana¬ 
lyzed in terms of multiple isotropic winds (Winters et al. 2003, 
2007). The interferometer data show an extended source of about 
15" (EWHM), and evidence for a ring structure in CO(2-l). 
However, when assuming multiple isotropic winds. Winters et 
al. (2007) did not manage to fully explain the line profiles ob¬ 
tained at different positions of the spectral maps. 

Composite CO line profiles have been observed for several 
AGB sources, including RS Cnc and X Her (e.g. Knapp et al. 
1998). Kahane &. Jura (1996) interpreted the narrow line com¬ 
ponent of X Her in terms of a spherically expanding wind and 
the broad line component in terms of a bipolar flow, probably 
observed at a viewing angle of about 15°. For RS Cnc, high spa¬ 
tial resolution observations can be interpreted assuming the same 
geometry but with a different inclination of about 40° with re¬ 
spect to the line of sight (Hoai et al. 2014). These results raise 
the question whether the bipolar outflow hypothesis might also 
apply to EP Aqr. 

The present study is therefore revisiting the analysis of 
Winters et al. (2007), assuming instead a bipolar outflow ap¬ 
proximately directed along the line of sight. We exploit the ex¬ 
perience gained in similar studies of other stars, RS Cnc (Hoai et 
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al. 2014, Nhung et al. 2015) and the Red Rectangle, apost-AGB 
source, which could be interpreted by applying a similar model 
(Tuan Anh et al. 2015). 

The paper is organised as follows: in Sect. |2]we present re¬ 
processed data on EP Aqr that will be used in this work. Sect. [2 
describes in an elementary way the main features implied by the 
assumed orientation of EP Aqr, with the star axis nearly paral¬ 
lel to the line of sight, mimicking spherical symmetry. Sect. |4] 
presents an analysis of the data using a simple bipolar outflow 
model that had been developed earlier for RS Cnc (Hoai et al. 
2014, Nhung et al. 2015); the results of the best fit to the obser¬ 
vations is used as a reference in the following sections. Sect. |5] 
presents a study of the CO(l-O) to CO(2-l) emission ratio, allow¬ 
ing for an evaluation of the temperature distribution nearly inde¬ 
pendent of the model adopted for the description of the gas den¬ 
sity and velocity. Sect. |6]uses the distribution of the gas velocity 
obtained in Sect. |4]to reconstruct directly in space the effective 
gas density, without making direct use of the parametrization of 
the flux of matter obtained in Sect. 01 and providing therefore 
a consistency check of the results of the model. Einally, Sect. |7] 
presents a detailed study of the spatial variation of the radial ve¬ 
locity of the narrow line component. It provides a sensitive check 
of the bipolar outflow hypothesis and of its small inclination an¬ 
gle with respect to the line of sight. 


2. Observations 

The observations analysed here combine interferometer data 
obtained using the Plateau-de-Bure Interferometer with short¬ 
spacing data obtained using the Pico Veleta 30-m telescope. 
A detailed description of the original observations is given in 
Winters et al. (2007). '2cO(2-l) and '^00(1-0) spectral data 
have been obtained with a spatial resolution of ~\" and ~2" 
respectively and a spectral resolution of 0.1 km s '. At such 
spatial resolution, the two spectral components present in the 
single-dish line spectra are seen as originating from the same 
region. The images are virtually circularly symmetric and dis¬ 
play smooth variations with velocity and projected distance from 
the star. In Winters et al. (2007) the observed spectra were in¬ 
terpreted by a succession of spherically symmetric and short¬ 
spaced mass-loss events. These authors noted the presence of in¬ 
homogeneities in the spectral maps indicating a possible dumpi¬ 
ness of the circumstellar envelope. 

We have reprocessed the original data, correcting for an ar¬ 
tifact introduced by a velocity shift of 0.52 km s ' in the CO(l- 
0) data and recentering the maps on the position of the star at 
epoch 2004.0, correcting for its proper motion as determined by 
van Leeuwen (2007). The reprocessed channel maps are shown 
in Eigs.[T]and|2 

3. Observing a star along its symmetry axis 

The high degree of azimuthal symmetry observed in the chan¬ 
nel maps of both CO(l-O) and CO(2-l) emission (see Eigs. [1] 
and |2] Winters et al. (2007)) implies that a bipolar outflow, if 
present, should have its axis nearly parallel to the line of sight. 
In such a configuration, some simple relations apply between the 
space coordinates and their projection on the sky plane. We re¬ 
view them briefly in order to ease further discussion. We use co¬ 
ordinates (Eig. [2l centered on the star with x along the star axis, 
parallel to the line of sight, y pointing east and z pointing north. 
We assume rotation symmetry about the x axis and symmetry 
about the equatorial (y, z) plane of the star. Defining p the gas 


'^CO emission from EP Aqr 



Right Ascension 


Fig. 1. Channel maps in the *^CO(J=1-0) line (smoothed to a 
width of Ikms^'). Contours are plotted at 5, 10, 20, 30, 40cr 
(1 cr = 14mJy/beam). The synthesized beam is indicated in the 
lower left. 



Right Ascension 


Fig. 2. Channel maps in the *^CO(J=2-l) line (smoothed to a 
width of 1 kms '). First contour is plotted at 10 cr, the follow¬ 
ing contours start at 20 cr and are plotted in 20 cr steps (1 cr = 
16mJy/beam). The synthesized beam is indicated in the lower 
left. 


density, T the gas temperature and V the gas velocity in the star 
rest frame, having components Vx, Vy, V^, the following relations 
apply at point (x,y,z): p and T are functions of r = (x^ +y^ +z^)^ 
and of \a\ where a = arctan{xlR) is the latitude in the star refer¬ 
ence system, and R - (y^ + z^)^ is the projection of r on the sky 
plane. 

The velocity V takes a simple form when expressed in terms 
of meridian coordinates, R = ycosw + zsincu, x and j] = 
-y sin « -H z cos a>, at being the star longitude; its components 
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Fig. 3. Definition of coordinates. The (y,z) coordinates are ob¬ 
tained by rotation of angle oj in the sky plane about the x axis 
(line of sight and star axis, cf. Sect. |3]l. In the sky plane, y is 
towards East, and z towards North. 

are then functions of r and |a|. In case of pure radial expansion, 
Vr^ - 0, Vt = xVjr and Vs = RVIr while in case of pure rotation 
Er = V^, = 0. 

In the case of pure radial expansion with a known radial de¬ 
pendence of the velocity, the flux density / measured in a pixel 
(y, z) at velocity Vx is associated with a well dehned position in 
space, obtained from 

V = rVxlV = + R")^-VxIV = RVxiV^ - V^xT^ (1) 

Then, if one knows V, one can calculate 14 and therefore its 
derivative dVxIdx at any point in space. Dehning an effective 
density p^g as p^g = fdVxIdx, such that the observed flux in a 
given pixel be J fdVx - f Pegdx, one is then able to calculate 
it at any point in space. The effective density, dehned by this 
relation, is the product of the actual gas density and a factor ac¬ 
counting for the population of the emitting state, the emission 
probability and correcting for the effect of absorption. 

We illustrate the above properties with a simple example, 
close to the reality of EP Aqr. We assume exact rotation invari¬ 
ance about the star axis (parallel to the line of sight), a purely 
radial wind with a velocity that only depends on a 

V -Vg + Vi sin^a (2) 

(To at the equator and Vo + Vi at the poles) and an effective 
density (dehned such that its integral over the line of sight mea¬ 
sures the observed Hux density) inversely proportional to r^. 
Erom Vx - V sin a and Eq. © we obtain sin^a + (Vo/Vi) sin a - 
(VxIVi) - 0. This reduced cubic equation is explicitly solvable 
with sin O' = -(q - s)^ - {q + s)^, s - (q^ + p p - 5 V 0 /V 1 
and q — -\VxlV\. Eor a given value of the Doppler velocity Vx 
one can therefore calculate a independent of R: in any pixel, Vx 
and a are related in the same way, which is illustrated in Eig.|4] 
(upper panel). Similarly, Vx and r/R - 1/ cos a are also related 
in a universal way, independent of R (Eig.lH middle panel). The 
resulting velocity spectra are displayed in Eig. |4] (lower panel) 
for R-1, 2, 3, 4 and 5 (all quantities are in arbitrary units ex¬ 
cept for the velocities that are in kms^' with Vo = 2 kms^' and 
Vi = 8 kms^*). They simply scale as 1/R, rehecting the Ijp- 
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Fig. 4. Relation between the Doppler velocity Vx and the star lat¬ 
itude a (cf. Eq. ©, upper panel), the ratio rjR — \ jcosa (mid¬ 
dle panel) for the simple star model described in the text. Lower 
panel: velocity spectra obtained for the same model at R=l, 2, 3, 
4 and 5 (running downward). 


dependence of the effective density. In such a simple model, the 
flux ratio between the CO(l-O) and CO(2-l) is a constant, inde¬ 
pendent of R and Vx- 


4. Comparison of the observations with a bipoiar 
outflow modei 

The spectral maps extracted from the reprocessed data (see 
Sect.© are shown in Eig.|5] The spectra are displayed in steps of 
1" in right ascension and declination. The synthesized beam is 
3.53"xl.84" (PA=28°) forCO(l-O), and 1.67"x0.94" (PA=29°) 
for CO(2-l), respectively. 

To describe EP Aqr, we use the model of Hoai et al. (2014), 
which assumes that the wind is purely radial, free of turbulence, 
and in local thermal equilibrium. Eor solving the radiative trans¬ 
fer, a ray-tracing method that takes absorption into account is 
used (Hoai 2015). Moreover, the wind is assumed to be station¬ 
ary and is supposed to have been in such a regime for long 
enough a time, such that the radial extension of the gas vol¬ 
ume is governed exclusively by the UV dissociation of the CO 
molecules by interstellar radiation (Mamon et al. 1988) and does 
not keep any trace of the star history. Eor EP Aqr the CO/H abun¬ 
dance ratio is taken as 2.5 10“^ (Knapp et al. 1998), a value rep¬ 
resentative for an M-type star. The temperature is parametrized 
by a power law of r,T - but is independent of the stel¬ 

lar latitude a. The spectrum of flux densities in each pixel is 
calculated by integration along the line of sight, the tempera¬ 
ture dependent contributions of emission and absorption being 
respectively added and subtracted at each step. The wind veloc¬ 
ity V and the flux of matter /m are smooth functions of | sin a\ 
and allowance is made for a radial velocity gradient described by 
parameters Ai at the poles and A 2 at the equator. Assuming the 
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Fig. 5. Spectral maps centered on the star of the CO observations (black) and the best-fit model (blue). The CO(l-O) map is shown 
in the upper panel, CO(2-l) in the lower panel. Steps in right ascension and declination are 1". 
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Table 1. Best fit parameters obtained for the CO(l-O) and CO(2- 
1) data. A distance of 114 pc is adopted. 


Parameter 

Best fit value 

e 

13° 


144° 

cr 

0.3 

Vi 

S.Okrns-' 

V2 

2.0 kms^' 

Imx 

1.26 10-’' Moyr-'sr-' 

/m2 

0.49 10“’' Moyr-'sr-' 

4, 

0.52 

^2 

0.38 

To {r = 1") 

116K 

n 

0.77 


wind to be stationary, the density p is then dehned at any point by 
the relation M{a) — A-nfMia) - Anr^p{a, r)V{a, r) — const{a). 

The bipolarity of the flow is parametrized as a function of 
sin a using Gaussian forms centered at the poles; 

G = ex/ 7 [-i(sina' - l)^/cr^] + exp[-^(sina + l)^/cr^], (3) 

where cr is a parameter to be adjusted. This function is used 
to dehne the sin a dependence of both the wind velocity and the 
flux of matter: 

V = yiG(l - + y2(l - (4) 

and 

fn - /mi G + /m 2 ■ (5) 

The density then results from mass conservation: 

p = r-\fM,G + /M,)/[yiG(l - A,e-'-l^-^") + y 2 (l - 

( 6 ) 

For r-independent wind velocities, Ai = A 2 = 0, the velocity 
and flux of matter are nearly Vi + V 2 and /m, + /mj at the poles 
and eV\ + V2 and s/m, + /mj at the equator where s is the small 
positive value taken by G at the equator. 

A small inclination 0 of the star axis with respect to the line 
of sight, with position angle / with respect to the north, is made 
allowance for. The relations quoted in Sect.[3]for 6 = 0 are mod¬ 
ified accordingly. More precisely, the position of the star frame 
with respect to the sky plane and line of sight is now defined 
by angles 0 and / and the polar coordinates in the star frame 
are the latitude a and the longitude cj defined in such a way as 
to conform with the definition given in Sect. [3] for 0-0. The 
star radial velocity as determined from the CO line profiles is 
-33.5 kms“* in good agreement with the value of -34 kms“* 
derived in Winters et al. (2003). The values taken by the param¬ 
eters are obtained by a standard minimization method. The 
best fit values of the parameters of the model to the joint CO(l- 
0) and CO(2-l) observations are listed in Table [1] and the quality 
of the fit is illustrated in Fig. |5] where modeled velocity spectra 
are superimposed over observed ones. 

This result shows that it is possible to describe the morphol¬ 
ogy and kinematics of the gas envelope of EP Aqr as the combi¬ 
nation of a slow isotropic wind and a bipolar outflow, the wind 
velocities at the equator and at the poles being ~ 2 kms * and 
~ 10 km s ' respectively and the inclination of the star axis with 
respect to the line of sight, 0, being ~13°. The flux of matter 
varies from 0.49 10“® Mq yr^' sr"' in the equatorial plane to 1.26 
10“* Mq yr^'sr^* in the polar directions. The total mass loss rate 



Fig. 6. CO(l-O) to CO(2-l) flux ratio (black) where each of the 
CO(l-O) and CO(2-l) fluxes has been averaged over 5 successive 
velocity bins and over the concentric rings defined in the text. 
The lower right panel is for all pixels having R < 5". The result 
of the best fit of the model described in Sect. |4] is shown in red 
and that of its modification described in Sect.|5]is shown in blue. 
The horizontal lines indicate the level (0.063) above which the 
data must be under the hypothesis of local thermal equilibrium. 

(integral of M(a) over a) is 1.8 10“’ Mq yr“*. The angular aper¬ 
ture of the bipolar outflow, measured by the parameter cr - 0.3 
in Eq. (O, corresponds to a = arcsin(0.7) ~ 45°. 

Leaving for Sect. Q a more detailed discussion of this result, 
we study in the next section the CO(l-O) to CO(2-l) flux ratio, 
which gives important information on the temperature distribu¬ 
tion with only minimal dependence on the details of the model. 

5. The CO(1-0) to CO(2-1) flux ratio 

Comparing the fluxes associated with CO(l-O) and CO(2-l) 
emission provides information on the gas temperature nearly 
independently from the gas density and velocity, which are 
common to the two lines. In the optically thin limit and as¬ 
suming local thermal equilibrium, the temperature ^(K) is ob¬ 
tained from the ratio Q of the CO(l-O) to CO(2-l) fluxes as 
T - ll.l/ln(g/0.063). However, the model presented in the 
preceding section does not reproduce well the observed flux ra¬ 
tio. In particular, we observe that Q is sometimes below the min¬ 
imal value 0.063 allowed by the local thermal equilibrium hy¬ 
pothesis (when T — 00 ). This is illustrated in Eig. | 6 ] where the 
data have been averaged over concentric rings limited by circles 
having R-V , 2", 3", 4" and 5" respectively. 

Accordingly, we allow for a small difference of normaliza¬ 
tion between the CO(l-O) and CO(2-l) fluxes by multiplying the 
former by a factor \ + p and the latter by a factor \ - p. This 
factor may account for a real difference in calibration as well as, 
in an ad hoc way, for any inadequacy of the model (assumption 
of local thermal equilibrium, absence of turbulence). Moreover, 
we see from Eig. that one needs to increase the ratio Q near 
the extremities of the velocity spectrum more than in the mid- 
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Fig. 7. Distribution of the gas temperature. Left; as a function of 
a at distances from the star r = 1" to 8" (from top to bottom) in 
steps of 1"; right; as a function of r at latitudes a = 0° (red), 30° 
(black), 60° (green) and 90° (blue). 


Table 2. Best fit parameters of the CO(2-l) to CO(l-O) ratio. 


Parameter 

FofK) 

a 

til 

111 

//(%) 

X^ldof 

Best fit values 

170 

0.67 

1.22 

-0.55 

9.0 

1.14 

Uncertainties 

20 

0.12 

0.06 

0.10 

3.0 


die, meaning near the poles at large values of r more than near 
the equator at small values of r. Accordingly, we parametrize the 
temperature in the form 

r = 7’o(sin^a + acos^a)(r/l")^” (7) 

with 

n - fii + n2 cos^a. (8) 

At r = 1", the temperature is To at the poles and aTo ^t the 
equator. The power index describing the radial dependence of 
the temperature is at the poles and «i + U 2 at the equator. 

The result is summarized in Table|2]and illustrated in Fig.|7] 
The ratio between the values taken by the temperature at the 
equator and at the poles, ar^"-, varies from 0.70 at r = 1" to 1.64 
at r = 5", crossing unity at r » 2". As r/R is much larger near 
the poles than at the equator (it is equal to l/cosa in the simple 
configuration studied in Sect. |3]l, the extremities of the velocity 
spectrum probe large r values, meaning low temperatures, as re¬ 
quired by the measured Q values; both the values taken by ^2 
and by a are determinant in increasing Q at the extremities of 
the velocity spectrum. For the same reason, the temperature in¬ 
crease obtained from the parametrization at small values of r and 
large values of a is not probed by any of the data, which makes 
its reliability uncertain. 

The values of the temperature that we find close to the central 
star (~ 1000 K at 0.1") are consistent with those obtained by 
Cami et al. (2000) using infrared CO 2 lines. 

The uncertainties listed in Table|2]are evaluated from the de¬ 
pendence of over the values of the parameters, taking proper 
account of the correlations between them. However, they have 
been scaled up by a common factor such that the uncertainty 
on Tq be 20 K. This somewhat arbitrary scaling is meant to 
cope with our lack of detailed control over systematic errors 
and the value of 20 K is evaluated from the robustness of the 
results against changes of different nature that have been made 
to the model finally adopted in the course of the study. While 
the quoted uncertainties give a good idea of the reliability of the 
result and of its sensitivity to the values taken by the param¬ 
eters, both the quality of the data (in particular the need for a 


//-correction) and the crudeness of the model do not allow for a 
more serious quantitative evaluation. 

We noted that the quality of the fits can be slightly improved 
by allowing for a temperature enhancement at mid-latitudes and 
small distances to the star; however, we were unable to establish 
the significance of this result with sufficient confidence. 

The following two conclusions can then be retained; 

- At large distances to the star, the temperature is significantly 
higher at the equator than at the poles, by typically 15 K at 
r = 5"; 

- The temperature decreases with distance as ~r at the 
poles and significantly less steeply at the equator, typically 


6. Evaluation of the effective density in the star 
meridian plane 

Under the hypotheses used to evaluate the model described in 
Sect, m one can associate to each data-triple (y, z, Vx) a point 
(r, a) in the meridian plane of the star. Indeed, .x is a known func¬ 
tion of r, a, and the star longitude ai, as is 14 = Vxj r. The former 
uses the values of 0 and ij/ and the latter the parametrization of 
y as a function of r and a that were obtained from the best fit. 
As one has three equations and three unknowns one obtains r, a 
and a> given Vx, y and z- The flux density / corresponding to a 
given data-triple can therefore be mapped as an effective density 
PefF = fdVx/dx in the meridian plane of the star, (i.e. each plane 
that contains the star’s polar axis) with coordinates r cos a and 
rsin a (see Fig.O. Its evaluation uses the best fit parametriza¬ 
tion of the wind velocity but not that of the flux of matter; it 
provides therefore a consistency check. 

The validity of the procedure has been checked by replac¬ 
ing the data by the results of the model and verifying that the 
reconstructed effective density is identical to that used as input. 
The region of large lal values is probed exclusively by the central 
pixels near the extremity of the spectrum and is therefore subject 
to large systematic uncertainties. Indeed, large \a\ values mean 
large |y;i| values, which in turn mean large r/R values, namely R 
much smaller than r. For this reason, on the basis of the result of 
the validity check, we restrict the space reconstruction to pixels 
having R > 1". In addition, in order to avoid dealing with too 
low fluxes, we require R not to exceed 10". The resulting maps 
of r^Peff are displayed in Fig. 0 for CO(l-O) and CO(2-l) sepa¬ 
rately and its averaged distributions as functions of r, a and a> in 

Fig. 121 

All distributions display significant deviations with respect 
to the distributions generated by the model, which, however, re¬ 
produces well the general trends. In principle, one should be 
able to include in the model the result of the space reconstruc¬ 
tion of peff and iterate the fitting procedure with a more realistic 
parametrization of the flux of matter and possibly a more sophis¬ 
ticated velocity model. However, the quality of the present data 
is insufficient to undertake such an ambitious programme. We 
shall be satisfied with simply taking note of the differences be¬ 
tween the model of Sect. |4] and the image revealed by the space 
reconstruction, as long as they can be considered as being both 
reliable and significant. 

The maps of Fig. [8] and the a distributions of Fig. |9] give 
evidence for a significant asymmetry between the northern and 
southern hemispheres of the star (defining north and south in the 
star coordinate system, north corresponding to a = 7r/2, away 
from Earth, with a projection towards North in the sky plane). 
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Fig. 8. Reconstructed maps of the effective density, multiplied 
by r^, in the star meridian plane under the assumption of a wind 
velocity having the form obtained from the best fit of the model 
of Sect.|4] The colour codes are such that the ratio between max¬ 
imum and minimum values of pr^ are the same for CO(l-O) 
(left) and CO(2-l) (right). Units are Jy beam^' kms^' arcsec. 
The rectangles show the regions selected for displaying the a> 
distributions shown in Fig. [10] (for CO(l-O), abscissa between 
0" and 3", ordinate between 5" and 9"; for CO(2-l), abscissa 
between 2.0" and 5.4", ordinate between 1.8" and 5.2"). 



Fig. 9. Distributions of the average effective density, multiplied 
by r^, reconstructed in the meridian plane of the star for a wind 
velocity having the form obtained from the best fit of the model 
of Sect. |4| Note that the yU-correction has not been applied, left; 
radial dependence of averaged over a and w; middle; lat¬ 
itude dependence of r^pes, averaged over r and w; right; longi¬ 
tude dependence of r^pefr, averaged over r and a. In each panel, 
CO(l-O) results are shown in red and CO(2-l) results in blue. 
The dashed curves show the results of the model, ignoring ab¬ 
sorption. 


Indeed, a modification of the model allowing for an asymme¬ 
try of the same form as used in Nhung et al. (2015) for RS 
Cnc, namely multiplying the flux of matter by a common fac¬ 
tor 1 H- y sin a, gives a 20% excess in the northern hemisphere 
of the star. Differentiating between the asymmetries of the two 
mass loss terms, /m, G and /mj reveals a strong correlation be¬ 
tween them, the asymmetry of the first term being more efficient 
at improving the quality of the fit. Restricting the asymmetry to 
the first term results in a 37% excess in the northern hemisphere. 
The effect of absorption does not exceed ~5% on average and 
cannot be invoked to explain such asymmetry. The fluctuations 



Fig. 10. Distributions of the effective density, multiplied by as 
a function of star longitude u for the annular regions delineated 
by the rectangles shown in Fig. [8] CO(l-O) data are shown in 
red, CO(2-l) data in blue. 


observed in the distributions of the effective density, as well as 
the differences observed between the CO(l-O) and CO(2-l) data, 
make it difficult to locate precisely the northern excess in space. 
Fig.dOjdisplays the u distributions of the effective density aver¬ 
aged over the r-a rectangles drawn in Fig.|8]around the maxima 
of the meridian maps; they do not show an obvious enhance¬ 
ment, the maximum observed in CO(2-l) corresponding in fact 
to a minimum in CO(l-O). This suggests that the northern ex¬ 
cess is distributed over the whole range of star longitudes rather 
than being confined to a particular region in space. This is con¬ 
firmed by the (y, z) sky maps of the measured fluxes (integrated 
over Doppler velocity and multiplied by R) that are displayed 
in Fig. [TT]and show an enhancement at radii exceeding ~3". 
This enhancement, which was already noted by Winters et al. 
(2007) in the case of the narrow velocity component for CO(2- 
1), is also visible in the broad velocity component. A better space 
resolution would be necessary to make more quantitative state¬ 
ments concerning this excess, its precise location, morphology 
and velocity distribution. The present data allow for a qualita¬ 
tive description, retaining as likely the presence of a northern 
enhancement at distances from the star exceeding 3" and dis¬ 
tributed more or less uniformly in star longitude. 


7. The mean Doppler velocity of the narrow line 
component 

The centroid of the narrow component of the observed spectra 
moves across the sky map by ~ + 0.5 km s“' on either side of the 
average velocity. We measure its mean velocity Av in each pixel 
after subtraction of the underlying broad component, interpo¬ 
lated from control regions spanning velocity intervals between 
5.4kms“' and 2.6kms“' with respect to the average velocity 
(either blue-shifted or red-shifted). The procedure is illustrated 
in Fig. [T2]for the spectra summed over the 37x37=1369 pixels 
of the sky map (9.25" x 9.25"). The resulting Av distributions, 
measured in km s ' with respect to the mean velocity measured 
for the whole maps, are shown in Fig. [13] (left) and their maps 
in Fig. [14] Both CO(l-O) and CO(2-l) data display very similar 
features, the narrow component being seen to be red-shifted in 
the north-west direction and blue-shifted in the south-east direc¬ 
tion, nearly parallel to the direction of the projected star axis. 
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y (arcsec) 


Fig. 11. Sky maps of Rf for CO(l-O) (upper panels) and CO(2-l) (lower panels) observations. Units are Jy beam ' km s ' arcsec. 
From left to right: all velocities, |14| < 2 km s ' and |14| > 2 km s'. The circles at /? = 3.5" correspond to the enhancement seen by 
Winters et al. (2007) in CO(2-l) data restricted to the narrow velocity component. The projection of the star axis on the sky plane 
(making an angle of 13° with the line of sight) and the axis from which the star longitude u) is measured (positive clockwise) are 
shown as black arrows. 



Fig. 12. Illustration of the procedure used to evaluate the mean 
Doppler velocity Av of the narrow component. The velocity 
spectra summed over the 1369 pixels of the map are shown in 
black for CO(l-O) (upper panel) and CO(2-l) (lower panel). The 
quadratic fit over the two control regions is shown in red and its 
interpolation below the narrow component in blue. The vertical 
blue lines show the mean Doppler velocities of the narrow com¬ 
ponent from which the interpolated broad component has been 
subtracted. They are used as reference for evaluating Av in each 
pixel separately. 


This suggests that the inclination of the star axis with respect 
to the line of sight, measured by the angle 0~13°, might be the 
cause of the effect. 

In order to assess quantitatively the validity of this interpre¬ 
tation, we define a band running from south-east to north-west as 
shown in Fig. [141 Calling f the sky coordinate running along the 


band, we display in Fig.[T3](right) the dependence of Av on ^ for 
pixels contained inside the band. It reveals a complex structure 
with successive bumps superimposed over a global increase, dis¬ 
playing remarkable similarities between CO(l-O) and CO(2-l), 
implying that it is likely to be real rather than an artifact of the 
analysis or instrumental. As a model, we use the best fit intro¬ 
duced in Sect. H) however with the temperature distribution and 
the adjustment of the relative CO(l-O) to CO(2-l) normalization 
that were introduced in Sect. [5] As illustrated in Figs. [TSl and [T6l 
the adopted model reproduces correctly the main features of the 
observed distributions: the width of the Av distribution and the 
amplitude of the Doppler shift when moving from south-east to 
north-west along the ^ axis. Again, as in the preceding section, 
observations are well reproduced by the model, not only quali¬ 
tatively but also quantitatively in their general features. This is 
only achieved by the hypothesis of a bipolar wind structure and 
is inconsistent with the assumption of a spherical outflow. 

8. Discussion 

Applying a bipolar outflow model to analyze the spatially re¬ 
solved CO(l-O) and CO(2-l) spectra from EP Aqr results in an 
excellent agreement of the model results with the observations. 
The peculiar geometry of the star orientation, with the symme¬ 
try axis nearly parallel to the line of sight, implies an approxi¬ 
mate circular symmetry when projected on the sky and gives the 
illusion of a spherically symmetric gas distribution. However, 
when correlated with the Doppler velocity distributions, it al¬ 
lows for a transparent interpretation in terms of the spatial ef¬ 
fective density, each data-triple being associated with a single 
point in space once the velocity of the radial wind is known. 
Under the hypotheses of pure radial expansion and of rotation 
invariance about the star axis, it is possible to reconstruct the 
space distribution of the effective density. As a first step, we have 
compared the observations with a simple bipolar outflow model 
that had been developed earlier to describe another AGB star. 
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Fig. 14. Sky maps of the mean Doppler velocity Av (kms“^) measured with respect to its value averaged over the whole map 
(shown as ’’reference” in Fig.fTSli for CO(l-O) (left) and CO(2-l) (right). The black lines limit the bands in which pixels are retained 
to evaluate the ^ dependence of Av (Fig. fOl right). 
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Fig. 16. Results of the model introduced in Sect. |4]for the maps displayed in Fig. [14] The red arrows indicate the projection of the 
star axis on the sky plane. 




Fig. 13. Left; distributions of the mean Doppler velocity Av 
(kms"') measured with respect to its value averaged over the 
whole map (shown as ’’reference” in Fig. [T2l i. The black curve 
shows a two-Gaussian common fit to the two distributions. 
Right: dependence of the projection of the mean Doppler ve¬ 
locity Av (kms *), averaged over pixels included in the bands 
shown in Fig. [141 on coordinate ^ measured from south-east to 
north-west. The black curve shows a polynomial common fit to 
the two distributions. In both panels the CO(l-O) data are shown 
in red and the CO(2-l) data in blue. 


RS Cnc. A satisfactory description has been obtained by adjust¬ 
ing the model parameters to best fit the CO(l-O) and CO(2-l) 
data together. Several parameters defining the morphology and 
the kinematics of the gas volume surrounding EP Aqr have been 
evaluated this way with good confidence. Among these are the 
orientation of the star axis, making an angle of ~13° with the 
line of sight and projected on the sky plane at ~144° from north 




Fig. 15. Results of the model introduced in Sect.[4|for the distri¬ 
butions displayed in Fig. [T3| The smooth curves are the results 
of the fits (respectively Gaussian and polynomial) made to the 
observations in Fig. [13] 


toward west. The velocity of the radial wind has been found to 
increase from ~2 km s“' at the star equator to ~10 km s“' at the 
poles. The flux of matter is also enhanced at the poles, however 
less than the radial wind, resulting in an effective density slightly 
enhanced in the equator region. This shows that an axi-symmetry 
such as that observed in EP Aqr (or RS Cnc, see figure 6 in Hoai 
et al. 2014) is more clearly revealed by velocity-resolved data. 
Eor instance infrared data obtained in the continuum, even at 
very high spatial resolution, could not reveal such an effect. Such 
a kind of morpho-kinematics could explain that asymmetries in 
AGB outflows are found preferentially through radio observa¬ 
tions obtained at high spectral resolution (e.g. Castro-Carrizo et 
al. 2010). However, this characteristic of the model is partly ar- 
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bitrary and its main justification is to give a good description of 
the observations. We do not claim that the form adopted in the 
proposed model is unique; on the contrary we are confident that 
other forms could have been chosen with similar success. In par¬ 
ticular, the hypothesis of pure radial expansion retained in the 
model should not be taken as evidence for the absence of rota¬ 
tion of the gas volume about the star axis. Indeed, such rotation, 
if it were present, would be virtually undetectable as the result¬ 
ing velocities, being nearly perpendicular to the line of sight, 
would not produce any significant Doppler shift. 

The assumptions made in formulating the model are always 
approximations, some time very crude, of the reality. In partic¬ 
ular, evidence has been obtained for significant departures from 
the symmetry with respect to the star equatorial plane assumed 
in the model, revealing an excess of emission at northern star 
latitudes over the whole longitudinal range at angular distances 
from the star exceeding ~3". This feature, at least part of it, had 
already been noted earlier (Winters et al. 2007) and thought to 
suggest the existence of successive mass loss events. In other 
words a succession of mass loss events can be mimicked by a 
latitudinal variation of the flux of matter. 

In general, while the proposed model has been very efficient 
at describing the main features and general trends of the flux 
densities, some significant deviations with respect to the obser¬ 
vations have also been found. The present data, in particular their 
spatial resolution, prevent a more detailed statement about the 
nature of these deviations. An important asset of the model is to 
illustrate the efficiency of the method used to reconstruct the spa¬ 
tial morphology and kinematics of the gas envelope. Once data 
of a significantly better resolution will be available, e.g. from 
ALMA, using the same methodology will reveal many of such 
details with unprecedented reliability and precision. 

The availability of measurements of the emission of two dif¬ 
ferent rotation lines of the same gas is extremely precious and 
has been exploited as much as possible, based on the current 
data. The flux ratio of two lines allows to evaluate the tempera¬ 
ture distribution of the gas in a region where the ratio between 
the populations and emission probabilities of the two states vary 
rapidly, which is the case for angular distances from the star 
probed by the present observations. Most of the details of the 
model are irrelevant to the value of the flux density ratio be¬ 
tween the two lines, which depends mostly on temperature under 
reasonable approximations (local thermal equilibrium, absence 
of large turbulence, etc.). Through this property, evidence was 
found for a significant enhancement of temperature at the star 
equator, therefore correlated with the orientation of the bipolar 
flow and at variance with the hypothesis of a pure radial depen¬ 
dence. 

A method allowing iterating the wind velocity and the effec¬ 
tive density in successive steps has been sketched. When data 
with significantly improved spatial resolution will be available, 
this will provide a very efficient analysis tool. The peculiar orien¬ 
tation of EP Aqr makes such an analysis much more transparent 
than in the general case where the star axis is neither parallel nor 
perpendicular to the sky plane. 

Furthermore, the good spectral resolution of the present ob¬ 
servations has made it possible to finely analyse the variations 
of the Doppler velocity of the narrow line component across the 
sky map, providing a sensitive test of the validity of the bipolar 
flow hypothesis and a quantitative check of the inclination of the 
star axis with respect to the line of sight. 

The present work leads us to prefer a bipolar, stationary, 
wind model over the spherical, variable, scenario that was pro¬ 
posed earlier. AGB stars have traditionally been assumed to be 


surrounded by almost spherical circumstellar shells. Deviations, 
sometimes spectacular, from spherical symmetry are observed 
in planetary nebulae (PNs). Such deviations are also observed in 
post-AGB sources and it is generally considered that they arise 
during the phase of intense mass loss, after the central stars have 
evolved away from the AGB. However, it has also been recog¬ 
nized that the geometry of circumstellar shells around AGB stars 
may sometimes show an axi-symmetry (e.g. X Her, Kahane & 
Jura (1996); RS Cnc, Hoai et al. (2014)). The origin of this axi- 
symmetry has not been clearly identified: magnetic field, stellar 
rotation, or presence of a close binary companion may play a 
role. Recently, Kervella et al. (2015) reported the detection of 
a companion to L2 Pup, another M-type AGB star, with a disk 
seen almost edge-on. The case of EP Aqr is particularly inter¬ 
esting because, with no detected Tc and a ratio ~ 10, it 

seems to be in a very early stage on the AGB. All these cases 
raise the interesting possibility that the deviations from spheric¬ 
ity observed in PNs and in post-AGB stars are rooted in their 
previous evolution on the AGB. 

9. Conclusion 

We have revisited the EP Aqr data obtained by Winters et al. 
(2007). We now interpret them in terms of an axi-symmetrical 
model similar to that developed for RS Cnc by Hoai et al. 
(2014). The expansion velocity varies smoothly with latitude, 
from ~2kms * in the equatorial plane to ~10 kms * along the 
polar axis, which for EP Aqr is almost parallel to the line of 
sight. The mass loss rate is estimated to ~1.8 10“^ Mq yr"'. The 
two stars look very similar, the differences in the observations 
being only an effect of the orientations of their polar axis with 
respect to the line of sight. There is evidence for a temperature 
enhancement near the star equator and a faster decrease along 
the polar axis than in the equatorial plane. 

Both RS Cnc and EP Aqr show peculiar CO line emission 
whose profiles consist of two components (Knapp et al. 1998, 
Winters et al. 2003). This may suggest that the other sources with 
such composite CO line profiles (e.g. X Her, SV Psc, etc.) might 
also be interpreted with a similar model, rather than with suc¬ 
cessive winds of different characteristics as has been proposed 
earlier. 
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